Search for: All records

We aim to automatize the identification of collective variables to simplify and speed up enhanced sampling simulations of conformational dynamics in biomolecules. We focus on anharmonic low-frequency vibrations that exhibit fluctuations on timescales faster than conformational transitions but describe a path of least resistance towards structural change. A key challenge is that harmonic approximations are ill-suited to characterize these vibrations, which are observed at far-infrared frequencies and are easily excited by thermal collisions at room temperature. Here, we approached this problem with a frequency-selective anharmonic (FRESEAN) mode analysis that does not rely on harmonic approximations and successfully isolates anharmonic low-frequency vibrations from short molecular dynamics simulation trajectories. We applied FRESEAN mode analysis to simulations of alanine dipeptide, a common test system for enhanced sampling simulation protocols, and compare the performance of isolated low-frequency vibrations to conventional user-defined collective variables (here backbone dihedral angles) in enhanced sampling simulations. The comparison shows that enhanced sampling along anharmonic low-frequency vibrations not only reproduces known conformational dynamics but can even further improve sampling of slow transitions compared to user-defined collective variables. Notably, free energy surfaces spanned by low-frequency anharmonic vibrational modes exhibit lower barriers associated with conformational transitions relative to representations in backbone dihedral space. We thus conclude that anharmonic low-frequency vibrations provide a promising path for highly effective and fully automated enhanced sampling simulations of conformational dynamics in biomolecules. 

Peptides that bind to inorganic materials can be used to functionalize surfaces, control crystallization, or assist ininterfacial self-assembly. In the past, inorganic-binding peptides have been found predominantly through peptide library screening. While this method has successfully identified peptides that bind to a variety of materials, an alternative design approach that can intelligently search for peptides and provide physical insight for peptide affinity would be desirable. In this work, we develop a computational, physics-based approach to design inorganic-binding peptides, focusing on peptides that bind to the common plastics polyethylene, polypropylene, polystyrene, and poly(ethylene terephthalate). The PepBD algorithm, a Monte Carlo method that samples peptide sequence and conformational space, was modified to include simulated annealing, relax hydration constraints, and an ensemble of conformations to initiate design. These modifications led to the discovery of peptides with significantly better scores compared to those obtained using the original PepBD. PepBD scores were found to improve with increasing van der Waals interactions, although strengthening the intermolecular van der Waals interactions comes at the cost of introducing unfavorable electrostatic interactions. The best designs are enriched in amino acids with bulky side chains and possess hydrophobic and hydrophilic patches whose location depends on the adsorbed conformation. Future work will evaluate the top peptide designs in molecular dynamics simulations and experiment, enabling their application in microplastic pollution remediation and plastic-based biosensors.